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Abstract 

We attempt to give a bird's eye view of the physical mechanisms 
leading to anomalous relaxation, and the relation of this phenomenon 
with anomalous diffusion and transport. Whereas in some cases these 
two notions are indeed deeply related, this needs not to be the case. 
We review several models for stretched exponential relaxation (dif- 
fusion in traps, broad distribution of relaxation times, two-step re- 
laxation) and insist on the physical interpretation to be given to 
the stretching exponent 0. We then discuss compressed exponentials 
which have been recently observed in a variety of systems, from soft 
glassy materials to granular packs. We describe a model where slow 
rearrangement events occur randomly in space and create long ranged 
elastics strains, leading to = 3/2. 

1 Introduction 



Normal diffusion usually refers to plain Brownian motion, and describes situations 
where the mean-square displacement of a particle (or the variance of the (log)-price 
of a financial asset) grows linearly with time. This property is in general only valid 
when time is sufficiently large compared to a microscopic correlation time. The 
diffusion law R 2 oc Dt is indeed tantamount to a complete decorrelation between 
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successive displacements. Anomalous diffusion, on the other hand, describes all 
other cases, where the normal diffusion law fails to describe the data, at least over 
some time interval [i<,i>] where the logarithmic slope 2u(t) = din R 2 /dint is 
observed to be either smaller than one (subdiffusion) or larger than one (superdif- 
fusion) pQ. Of course, the effect is only interesting if the time interval [t < ,t > ] is 
wide enough to delimit a genuinely anomalous regime, and not a trivial crossover. 
One mechanism leading to anomalous diffusion is the existence of persistent (or 
anti-persistent) correlations up to time interval t>. Positive correlations of dis- 
placements lead to superdiffusion, whereas negative correlations lead to subdif- 
fusion. For example, diffusion of particles in hydrodynamical flows often lead to 
superdiffusion since particles are coherently convected along streamlines over some 
length scale p~j. Conversely, particles in a potential field tends to be trapped in 
minima, where they execute numerous back and forth oscillations, and therefore 
anti-correlated motion, before being able to escape pQ. Another important mecha- 
nism generating superdiffusion is the existence of jumps with a broad distribution 
of possible sizes, leading to "Levy flights" [2]. 

Normal diffusion is also characterized by a Gaussian diffusion profile, i.e., the 
distribution of the displacement R between to and to + t is given by: 



P(R, t) = , 1 exp 



4Dt 



(1) 



The relaxation of a density fluctuation of wavevector q, created at time to, is easily 
computed to be exponential^ 

C(q,t) = ^ (j2 ex P (»j(<0 + *) ~ rj(t ))]J = eM~Dq 2 t), (2) 

with relaxation time \/Dq 2 equal to the diffusion time over the scale 1/q. There- 
fore, normal diffusion is associated to what is usually considered to be "normal" 
- i.e. exponential - relaxation. The aim of this short review is to discuss possi- 
ble forms of anomalous relaxation in complex systems, in particular the relation 
between anomalous diffusion and non exponential relaxation. 



-'^Here we imagine that particles do not interact, and that there is no difference between 
the individual and collective diffusion constant D. More generally, there is a difference 
between the relaxation of tracer particles (of small concentration) and of the density of 
interacting host particles. 
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2 Stretched exponential relaxations 



Stretched exponential relaxation refers to cases where the correlation function 
decays as In C(t) oc — 1@ with (3 < 1. A simple model where this behaviour is 
exactly observed is the case where subdiffusion is induced by an anti-persistent 
Gaussian process (also called fractional Brownian motion [3]), in which case the 
diffusion profile is given by: 



P(R,t) 



exp 



R 2 



(3<l. 



(3) 



The density relaxation is then trivially given by C(q, t) = exp(— Dpq 2 t@). However, 
it is hard to think of a physical system for which this model applies directly. A 
more generic model for subdiffusion is anomalous trapping, for which, as we discuss 
now, the correlation indeed decays as a stretched exponential but only at small 
times. 



2.1 Trapping induced subdiffusion 

We therefore consider the trap model where a particle hops on a regular lattice in 
d dimensions; each node of the lattice is an energy valley of random depth E which 
traps the particle for a certain random time, before an activated event occurs and 
allows the particle to jump to a neighbouring site. One should distinguish the 
annealed case where the depth E is drawn anew at every jump of the particle 
from the quenched case where the depth E of each site is drawn once and for all 
and does not evolve with time. However, in d > 2, the particle typically visits 
each trap a finite number of times and the difference between the two models is 
immaterial at long times. We will therefore focus on the first case, which is much 
easier to solve analytically, and return to the quenched case in d = 1 at the end of 
this section. 

Conditioned to being in a trap of depth E and at temperature T, the exit time 
r is distributed according to: 

9 E (r) = T e-? E e"( r ° e ~" B >, (4) 

where Tq is a microscopic frequency and = 1/T (we set ks = 1 throughout the 
paper). If p{E) denotes the a priori probability distribution of the depth E of 
the traps, the particle performs a random walk among traps with a distribution of 
trapping times 

poo 

*(t) = / dE p{E) ^> e {t) (5) 
Jo 
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We will consider below the case where p(E) = 1/T g exp(—E/T g ) (exponential 
model), for which a true dynamical phase transition appears: diffusion is asymp- 
totically normal for T > T g , but subdiffusion set in for T < T g , with an anomalous 
diffusion exponent v (defined in the introduction) given by v = T/2T g < 1/2 . 
The exponential model is such that the upper limit t > of the time domain where 
subdiffusion is observed is infinite. Other cases could be considered as well, such 
as a Gaussian distribution of depths, but in that case t> would be finite and a 
crossover to normal diffusion would eventually set in. However, the following con- 
clusions will still hold in an intermediate time domain, which can become very 
large when T — > (see the detailed discussion in [I] ) . 

Interestingly, the subdiffusive phase T < T g of the exponential model is also 
non stationary, in the sense that all dynamical properties depend on the age of 
the system [6j 01 [5]. More precisely, if the initial (to = 0) position of the particle 
is chosen randomly, the two-time correlation function 



C(q,t + t w 

1 b w ) 



,iq- (f(t+t w )-r(t w )) 



(6) 



strongly depends on t w for T < T g , but becomes time translation invariant and 
given by Eq. ([2]) for T > T g , at least when Tot w 3> 1. The aging properties of the 
low temperature phase are strongly reminiscent of the phenomenology of glasses 
and spin-glasses below their glass transition [7], hence justifying the notation T g 
for the critical temperature of the model. 

The exact computation of C(q,t + t w ,t w ) was presented in [4]. After a se- 
ries of manipulations, one obtains the following result for the Laplace transform 
C(q,X,t w ) = f™dt e~ xt C(q,t + t w ,t w ): 



C(q,X,t w ) =n(\,t w ) 
where we have introduced: 



1 ~ Cd(q) 

l - *(A)Cd(g) 



+ 



1 1-*(A) 
A (CM)- 1 ~ *(A) 



(7) 



Cd(q) = j cos (^ c 



K=l 



and the following Laplace transforms: 



(8) 



*(A) 



dr e~ Xr *(t); U(X,t v 



dt e~ xt U(t + t w ,t h 



0) 



where H(t + t w ,t w ) is the probability not to have jumped between t w and t w + t. 
The above formulas are general and do not depend on the specific choice of p(E). 
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In the case of an exponential density of trap depths, and when fj, = T/T g < 1, the 
correlation function II is known and its Laplace transform reads [H d] : 



1 r° ■ -(Xt w )u f 1 , n _ 



IL(\,t w )~— / due~^J u - Z dv{l-vf-'v-^ (10) 

t w JO 7T 



l + u 



with the asymptotic behaviors 



1 1 

1 



r( M ) (A^)!-". ' (l2) 

Interestingly, the correlation function IT depends on the ratio t/t w when fi < 1: 
the characteristic time needed to leave a trap is therefore proportional to the age 
of the system. For /i > 1, on the other hand, one would find that II only depends 
on Tot and not on t w - the problem is time translation invariant in that case, as 
usual for equilibrium dynamics. 

Coming back to Eq. ((TJ), it is convenient to introduce the subdiffusion time t q 
defined as: 

(r„«,r = ^ d3) 

corresponding to the typical time needed by the particle to spread over a region 
of size 1/q. There are therefore three time scales in (|7|) : A -1 ,t g ,t w . We are 
interested in the region where all three are much larger than the microscopic time 
scale Tq 1 , where the results become universal in the sense that they only depend 
on the exponential form of p(E) at large E, and not on the details of p(E) for finite 
E. Still we have to distinguish various time regimes in ([7]). Let us first consider 
young systems, for which t w <C t q ,t. In this case, we find that at small times 
t <C t q , the correlation decays as a stretched exponential with exponent fx < lH 

, s sinc[7rul / t\ fM . t . 

]nC(q,t + t w ,t w )~- r ^ l + ^ (-) for^<t<i ? , (14) 

whereas at large times t^>t w , the correlation decays as a power-law (much slower 
than a stretched exponential): 

C(q,t + t w ,t w )~T(l+iJ,)(^) fortyj^tq^t (15) 



2 We define in the following sincrr = sinx/x. 
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In the aging regime, on the other hand, we find at small times: 

lnC(q,t + t w ,t w ) ~ — sinc[7r(l — //)] I — J for t q <C t <C t w , (16) 

, r,,- x shield i . , . 

\nC{q,t + t w ,t w ) ~ 1 ^ J — t fori<i to and t <g (17) 

and finally, at long times: 

C(q, t + t w ,t w ) ~ sinc[7r/i] for t q <^t w <^t (18) 

There are four interesting points to notice: 

• For qa very large, such that t q <C t, t w , we find the same asymptotic behaviors 
as for Il(i + t w ,t w ). Physically, this means that as soon as the particle 
has jumped once, the rapidly oscillating correlation function averages to 
zero. Hence, only the particles which have not yet moved contribute to the 
correlation. 

• There are two regimes where the correlation function behaves similarly to 
a stretched exponential, but only at small times, when t q <C t <C t w or 
t w <C t <C tq. The exponent (3 of this stretched exponential is however 
different in both cases: it is equal to /3 = fj, when t w <C t <C t q , and equal to 
f3 = 1 — [i in the other case. 

• In the regime t -C t w ,t q , we find an interesting sub-aging behaviour, where 
C(q,t + t w ,t w ) is a function of t/t]^^, see the discussions in |9l [8]. 

• The frequency dependent susceptibility defined by: 

POO _ 

X(q,v,t w ) = l + iu dt e tujt C(q,t w + t,t w ) = l + im C(q,—iu,t w ), (19) 



behaves in a Cole-Cole fashion for young systems: 

^■"■ f - =o '" 1+( -W (20) 

for ijj <C Tq, qa <C 1. In the aging regime u>t w 3> 1, the behaviour of % is 
given by: 

x(g> ' g " r(,)(-L m )^ (21) 

A similar expression was obtained in the context of spin-glasses in [6], and 
discussed further in the context of experimental data in [7]. Note that the 
noise spectrum S(uj) of the system is related to the susceptibility through a 
fluctuation-dissipation relation S(uj) oc \"{oj)/oj. 



6 



Two further remarks on this spatial trap model: a) the marginal case T = T g , 
H = 1 is quite interesting, since it exhibits exact dynamical ultrametricity |10| 
and weakly aging 1/f noise, see Eq. ([2T]) above and [6j[7]; b) the quenched case 
in d = 1 reveals a number of subtle peculiarities, discussed in [8]. In particular, 
the diffusion exponent v is now given by v = [J>/(1 + //)■ Using the results of [8], 
one finds that in the case of young systems, the result for the correlation function 
reads: 

\nC(q,t + t w ,t w ) ~ - + " ; C(q,t + t w ,t w ) ~ (j)", (22) 

for f <C t g and t 3> t g , respectively, with now t q given by: ga(rot g ) l/ = 1. The 
above shape is different from the annealed case, but the initial decay is still of the 
stretched exponential type. The asymptotic behaviour of C(q,t + t w ,t w ) in the 
aging regime can also be inferred from the results of [llj . We also refer to refs. 
[12j for more work on the trap model and some application to glassy dynamics. 



2.2 Two-step relaxation 

In the above model, there is no intra-node dynamics: traps are assumed to have no 
spatial extension. This is clearly unrealistic and the short-time, high-q behaviour 
of the correlation function should be sensitive to the intra-trap dynamics. As a 
simple model for this, we consider harmonic potential wells again organized at 
the node of a regular lattice. Each particle oscillates at the bottom of the well, 
and with a small probability per unit time proportional to A = Toexp(—E/T) 
the particle hops to its nearest neighbour site, where it lands at random with the 
equilibrium probability inside the harmonic well. We neglect for the moment the 
fluctuations of E, which is justified when T T g . The position of a particle at 
time t will be written as: 

f(t) = R(t)+5(t), (23) 

where R(t) labels the lattice site to which the particle 'belongs' at time t, and 5(t) 
is the position of the particle within the well, the center of the well being defined 
by S(t) = 0. The probability for the particle not jumping site between to and to + t 
is 

n(t) = exp(-At). (24) 

If the particle has not jumped, then <5(to + t) — 5 (to) is a Gaussian variable of 
variance given by: 

((5(t ) - 5(t + t)) 2 ) = 2A 2 [1 - exp(- 7 t)], (25) 
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where Ao is the width of the explored region of the well, and 7 the inverse relaxation 
time in the well. If on the other hand at least one jump has occurred between to 
and to + t, one has: 

((5(t )-5(to + t)) 2 ) = 2Al (26) 
Our aim is again to compute the correlation function: 

C(q,t) = ( e W*o)-r-T*o+*))}. (27) 

This quantity can be easily computed by separating the no jump situations from 
the situations where at least one jump has taken place. One finds: 

C(q, t) = n(t) e -9 2A o[ 1 - ex P(-^)] + (1 _ U(t))e- q2A o ^HR{t )-R{t +t)))j^ ( 28 ) 

where the subscript J indicates that at least one jump has taken place. Clearly, 
this last average is related to the unconditional average by: 

/ e *9-(-R(*o)--R(*o+t))\ = jjm + (1 — n(t))(e**"^* ^ - ^* 0+ *^)j. (29) 

Finally, the unconditional average is a classic result of random walk theory on 
lattices: 

( e i<HR(t )-R(to+t))} = exp [At(&(g) - 1)] . (30) 
In the limit qa -C 1, the final result reads (see [13] for similar calculations): 

C(q, t) = U(t) [ e -</ 2 Ag[l-cxp(- 7 t)] _ e -^ A §j + e - 9 2 Ag e -AM^ (31) 

Note that C(q,t = 0) = C{q = 0, t) = 1, as it should. Suppose At <C 1 (i.e. the 
probability of a jump is very small, II(i) ~ 1) and jt S> 1, corresponding to full 
equilibration within the initial well. In this limit, one then finds: 

C(q,t)^f[l + q 2 A 2 eM-7t)}, (32) 

i.e. an exponential convergence towards a plateau value / = e~ q2A o. If on the 
other hand At <C 1, many jumps have occurred (II(t) ~ 0) and: 

C(q,t) « fe- A{qa)2t . (33) 

The parameter / measures the relative weight of the two relaxation processes 
(intra-trap and inter-trap); and has the following interpretation: if one plots the 
above relaxation function in a log-lin representation, one observes two-step re- 
laxation, with a quasi-plateau of the correlation function, of height /. The fast, 
intra-trap process does lead to a complete decorrelation - on the contrary, for time 
scales much smaller than 7 _1 , the system appears to be non ergodic (it keeps a 
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non zero correlation with its initial state); the strength of this non ergodicity is 
measured by /. The quantity / is called the non-ergodicity, or Edwards- Anderson 
parameter, in the context of glassy dynamics. The above simple model also allows 
to discuss in very simple terms the non-Gaussian parameter a often discussed in 
this context. A purely Gaussian probability density is such that lnC(q,t) oc q 2 . 
The kurtosis a is in fact related to the coefficient of the q term of the expansion 
of In C(q, t) in the vicinity of q = 0. One finds that a(t) ~ [1 — H(t)]/t 2 for small t 
and a(t) ~ H(t)/t 2 at large t. In the naive model above, 1 — II(i) = At at small t, 
leading to a diverging kurtosis at small times. This is however an artefact of the 
instant 'jump' description of the way a particle exits from a trap. A more realistic 
model should account for the fact that the time needed to exit a trap cannot be 
infinitely small, leading to a much smaller value of 1 — II (t) at small t and there- 
fore a vanishing a(t) in that limit. Therefore, we expect a(t) to peak around the 
typical time needed to exit a trap, before decaying back to zero at larger times, as 
often seen in experimental or numerical data (see e.g. |14j). 

When the fluctuations of the depths E are reinstalled, in particular in the 
glassy phase T < T g , one finds that the above short time behaviour of C(q,t) is 
unaffected, while the long time behaviour is still given by the expression of the 
previous section, up to a factor / which describes the initial fast fall-off of the 
correlation function. 

2.3 Superposition of relaxation times and stretched 
exponentials 

In the models considered above, the stretching of the relaxation is associated to 
(sub-)diffusive entities. This needs not be the case; in some cases there is relaxation 
but no transport. Suppose for example that the physics is governed by two- level 
systems, rotors or dipoles which do not wander in space at all, but feel a ran- 
dom field which makes the local relaxation time r random as well. More precisely, 
suppose that there is a whole spectrum of independent, exponential relaxation pro- 
cesses, with a density of relaxation times p(r), and that each elementary relaxation 
process contributes equally to the correlation function, so that: 



Jo 

Suppose further that the resulting correlation function C(t) can be well fitted by 
a stretched exponential: 



what can we infer about the shape of the density of relaxation times p{r)7 After 
a very small time, the correlation function has dropped from its initial value of 1 




(34) 



C(i)«exp-(7i) 



(35) 
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by an amount AC ~ J p(r)dr: all relaxation times larger than t have not had 
time to relax yet, whereas all relaxation times smaller than t have already relaxed. 
Equating this result with the initial decay of the stretched exponential, one finds: 

f P {T)dT « ( 7 tf — » P (t) a T ^ 7V" 1 , (36) 
Jo 

showing an (integrable) divergence of the density of small relaxation times. 

Therefore, although stretched exponential relaxation is usually associated with 
slow relaxation, the above result implies a profusion of short time scales in the 
system. However, since the long time behaviour of the relaxation is much slower 
than an exponential, there must also be abundant long time scales. How abundant 
can be guessed by a saddle point calculation. Suppose that p(r) decays itself, for 
large r as exp(— Bt@ ). This implies, for large t, a decay of C(t) given by: 

up to sub-leading power-law corrections. For example, a significantly stretched 
relaxation = 0.5 corresponds to /3' = 1, i.e., an exponential decay of long relax- 
ation times [15J. Therefore, only a very limited density of very large relaxation 
times is enough to transform a pure exponential relaxation into a stretched one. 
The pure exponential case corresponds formally to f3' = oo, i.e. no long times at 
all. 

More generally, any non trivial distribution of relaxation times leads to a re- 
laxation function that is faster than exponential on short times and slower than 
exponential on large times, in the following sense: define the average relaxation 
time as (r) = C(t)dt = f °° tp(t)oIt; the reference exponential relaxation is 
taken to decay over this average time. Then, the initial slope of the relaxation is 
given by: 

A(~i roo 111 1 

~^7= -P^ dT = <-> ^ 7T' (Provided (-) < +cx>) (38) 
at Jo t t (r) r 

where the equality is reached only if there is a single relaxation time scale, p(r) = 
5(t— 7~ 1 ). Conversely, for large time scales, one can easily show that the relaxation 
can only be slower than exp(— t/(r)). 

An important case, which superficially corresponds to the trap model discussed 
above, is: 

Pir) = ^~Q{r - rol (39) 

i.e. a power-law density of relaxation times. The cut-off for small r (that could be 
chosen to be less abrupt) and the condition p > are needed for the distribution 
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to be normalisable. Because of the cut-off, (-) < +00; the initial slope of the C(t) 
is then finite and the short time decay is regular, at variance with the stretched 
exponential case described above. However, the long time decay is in that case 
much slower than stretched exponential: 



The above result is valid for all values of /i, in particular when [i < 1, which corre- 
sponds to an infinite average relaxation time. At this stage, the reader might feel 
completely nonplussed: why all the aging properties that we found within the trap 
model when /j, < 1 seem to have disappeared? The answer is that the two models 
are in fact fundamentally different. In the trap model, particles progressively dis- 
cover their environment and, as they explore space further, they fall in deeper and 
deeper traps which they never encountered before. When T < T g , the deepest trap 
encountered at time t w is so much deeper than all the others that it dominates 
the dynamics. In the second model, on the other hand, there is a relaxing entity 
for each and every trap, and therefore the whole distribution of relaxation time is 
probed right away. There is no aging whatsoever in that case. 

It should be noted that a stretched exponential fit to non exponential relaxation 
is often an acceptable fit of empirical (or numerical) data. For example, when r oc 
exp(E/T) with a Gaussian distribution of barrier heights, the relaxation obtained 
from Eq. f|34[) can still be approximatively fitted by a stretched exponential, with, 
as a rough rule of thumb [4] : 



where Eq is the width of barrier distribution. When a stretched exponential fit is 
attempted, it is important to specify whether it is supposed to be faithful to the 
initial decay of the relaxation t < 7 , and hence to probe the profusion of high 
frequencies in the system, or if it is rather the long time (t 3> 7 _1 ) aspect of the 
phenomenon which is targeted, which is often the regime where exact analytical 
results are available but where the correlation is so small — In C(t) 3> 1 that it 
nearly impossible to measure. To illustrate this point, is interesting to discuss 
the problem of survival of a Brownian particle in a random environment, where 
infinitely deep traps are randomly scattered in a d dimensional space. For this 
problem, a rigorous result, due to Donsker and Varadhan |16j . asserts that the 
total survival probability, averaged over the starting point of the Brownian particle, 
decays asymptotically as a stretched exponential, with j3 = d/(d + 2). This result 
can be simply understood using the above relation, Eq. ([37|) . between j3' and (3: in 




(40) 




T 



(41) 
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the model at hand, long times are associated with exceptionaly large regions free of 
traps where the particle can survive much longer than on average. These occur with 
probability Inp oc — L d , and the associated survival time is r(L) ~ L 2 /D, from 
which one immediately concludes that /?' = d/2. Although this result is rigorous, 
its domain of validity is, in d > 2, so far into the asymptotic regime that it is of 
little interest, at least for numerical or experimental purposes, because the number 
of surviving particles in this regime is exponentially small. For intermediate times, 
one finds that the survival probability can be fitted by a stretched exponential, 
but with an effective value of (3 7^ d/(d + 2). We should also mention that this 
survival model is in fact very subtle: although the average survival probability is 
given by the above Donsker-Varadhan result, the survival of a given particle with 
fixed initial positiion decays as — lnC(i) ~ t/lnt! We refer to the insightful papers 
of Ben Arous et al. |17l [T8] for a comprehensive account of this result. 



3 Models of compressed exponentials 

Again, a trivial model leading to both super diffusion and compressed exponential 
relaxation (i.e. /3 > 1) is when the noise driving the particle is a long-ranged 
persistent Gaussian process, in which case the diffusion profile is given by: 



P(R, t) = — ^=^^- exp 



R 2 



(3>1. (42) 



and therefore C(q,t) = ex.p(—D/3q 2 t^), which is faster than exponential when > 
1. As we noticed for anti-persistent Brownian motion, physical situations where 
this model applies directly are scarce. In fact, natural models for superdiffusion 
do not lead to compressed exponentials, as we now discuss. We then turn to the 
necessary ingredients needed to explain why compressed exponential relaxations 
have recently been observed in a variety of soft glassy materials - where one would 
have rather expected, because of their glassiness, stretched exponentials! On the 
other hand, exactly as stretched exponential relaxation means super-abundance 
of short relaxation times, compressed exponential relaxation requires that no fast 
relaxation channels exist. 



3.1 Superdiffusion and Levy flights 

The best known model of superdiffusion is a random walk with a distribution of 
jump size t which decays for large i as: 

m - |S (43) 



12 



with /x < 2 such that the variance of the distribution diverges. In this case, one 
knows that for large times, the distribution of the total displacement of the particle 
is given by a Levy stable law of index /xH 

P(R,t) oc ^L M (£); i/ = ± (44) 

where we have been sloppy with prefactors and constants. The crucial point is 
that, although has no simple expression in general (except for \i = 2 where it 
reduces to the Gaussian and fj, = 1 where it is the Cauchy distribution) , its Fourier 
transform is extremely simple. In fact, the relaxation of a density perturbation is 
an unspectacular exponential: 

C(q,t) =exp(-DqH). (45) 

Although diffusion is anomalous, relaxation is normal! The only anomalous feature 
is the scaling of the wavevector q, raised to the power fi rather than q 2 for normal 
diffusion. Exactly this behaviour was observed for tracer particles in a system 
of polydisperse, elongated micelles |19j . The very strong contrast between the 
diffusion constant of small micelles and long micelles was argued to be responsible 
for the broad distribution of jump sizes, and in turn for the Levy flight motion 
clearly observed experimentally |19j . 

Another well studied case of superdiffusion is random advection in long range 
correlated hydrodynamical flows. Consider the following Langevin equation for a 
particle in a d = 3 random flow: 

dx -* 

— = - V (x)+f}(x,t), (46) 

where fj is the usual Langevin noise and V is a divergence free random convection 
field, with long-ranged correlations: 

(ViWVjm ^ %,\x-y\- a , a<2, (47) 

\x— y\— >oo 

where T is a tensor ensuring that V is indeed an incompressible flow. This model, 
with a = —2/3 can be seen as a very rough model of tracer dispersion in a turbulent 
flow with Kolmogorov scaling- in order to be more realistic one should allow the 
flow field to become time dependent, with a correlation that includes a function of 
(t— t')/\x— y\ 2 / 3 . For the above static turbulence model, one can show that diffusion 
is anomalous with a diffusion exponent exactly given by v = 2/(2 + a) > 1/2 (see 

3 For simplicity, we restrict here to one dimensional Levy flights. 



13 



PQ), leading to the Richardson scaling v = 3/2 for a = —2/3. The shape of 
diffusion profile at large times is not exactly known. However, since the simplest 
self-consistent re-summation scheme of the perturbation theory for this problem 
(called Mode-Coupling or Direct Interaction Approximation (DIA) in different 
contexts [20J) gives the exact diffusion exponent u, one can hope that it also leads 
to an accurate shape of the diffusion profile. More precisely, the DIA suggests 
that the Fourier-Laplace transform of P(R, t) reads, in the long wavelength, low 
frequency limit pQ: 

C(£A)ra 1 —r^, (48) 

which is exactly the Fourier-Laplace transform of a Levy distribution! The in- 
teresting consequence is that diffusion in a long-range correlated hydro dynamical 
flow is equivalent, at long times and in the scaling regime R ~ t u , to a Levy flight 
(provided of course the DIA is trustworthy). But is also means that this model 
cannot generate compressed exponentials! 



3.2 Rearrangement-induced stress fields in elastic me- 
dia 

Still, one of the major experimental surprise of the last decade is the discovery 
of compressed exponential relaxation in a bevy of different soft matter materials, 
ranging from polystyrene micro-sphere gels, diblock copolymers, laponite, etc. [21] . 
These systems have very slow dynamics and rheological properties typical of glasses 
or jammed granular media, which would have suggested that relaxation functions 
would be found, as in glasses, to be stretched exponentials. More precisely, one 
finds experimentally that the correlation function decays as: 

C(q,t) = exp[-A(qt)% (49) 

with /3 ~ 3/2 > 1. A notable feature of the above result is the scaling between 
space and time q~ l ~ t, characteristic of ballistic effects. It was suggested early on 
[22] that such a result might be related to the elasticity of these materials which 
mediate the long-ranged stress fields generated when a local rearrangement occurs. 
In this section, we show how a more detailed model of the deformation induced by 
these local events can give rise to anomalous q and t behaviour of the structure 
function [23J. 



The model 

Following [22], we assume that the dominant mechanism is the random appear- 
ance of localized rearrangements. For example, in the system studied in [22 J, the 
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micro-particles forming the gel attract each other rather strongly, the gel tends to 
restructure locally so as to create dense packings of particles. Since the collapsing 
particles belong to a gel network, their motion will induce a certain strain field 
around them; other particles therefore move and dynamical light scattering probes 
this motion. Similarly, in other soft elastic media, any local rearrangement will 
induce a strain that propagates into the system. A rearrangement of particles will 
result in a force dipole of intensity Pq in a certain direction n. When the dipole is 
formed at point ro, the elastic strain field u at point f can be computed assuming 
a simple central force elasticity: 

KAu = -f(f), (50) 

where K is a compression modulus and f(r) is a local force dipole of the form: 

f{r) = f [e ■ V5(f - r )] n, (51) 

where we will assume that e = eft, i.e, the mean displacement of the particles 
creates a force in the same direction. This dipolar force can be simply expressed 
by its Fourier transform: 

f(k) = iPo (k.n) n exp(-ik.fo) (52) 

with Pq = /oeQ The solution of equation (|50|) is of course: 

■+/-A P ° (r-ro)-n^ 

4-7TA \r — ro\ 

The r -2 dependence of the strain field has an immediate consequence which will be 
of importance in the following: if there is a finite density of force dipoles randomly 
scattered in space, the probability that the stress at a given point has an amplitude 
u decays for large u as ii _5//2 , which has a diverging variance. (This divergence is 
however cut-off if the finite size of the dipoles is taken into account). This property 
of the distribution of displacements and stresses will be responsible for the unusual 
^-dependence of the correlation function. 

Now, let us assume that the rearrangement events are not instantaneous but 
take place slowly, over a certain time scale 9. This will turn out to be crucial in 
the following. A given event j starts at time tj and is completed at time tj + 8; 
the dipole intensity P(t) at time t is a certain function of (t — tj)/6, which can be 
taken to be approximately linear: Pj(t) = Po(t—tj)/8 before saturating at Pq when 



4 Note that a more refined model with shear modulus could be also be considered, but 
would only change some numerical factors in the following calculations. 
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the collapse is completed. The dynamics of the individual particles is presumably 
dominated by viscous friction, therefore we write the following equation of motion 
for the strain field u: 



du(f, t) 



KAu + Y,fj(r,t) + ff(f,t) 



(54) 



7 



at 



3 



where the Fourier transform of the dipolar force is 



fj(k,t) = iPj(t) (k.fij)fij exp(-ik.fj) 



(55) 



and 7 is a friction coefficient, and ff is the thermal random force due to the viscous 
bath, uncorrelated in time and in direction. We will finally assume in the following 
that these events occur randomly in space and time, with a certain rate p per unit 
volume and unit time, and with a random orientation of the force dipole n. The 
quantity K/j is a diffusion constant that we will call D. Equation (|54|) defines 
the model that we want to study and from which we will compute the correlation 
function (or dynamical structure factor) C(q,t), defined as: 



where the brackets refer to a spatial average over r or, equivalently, over the 
random location and time of the rearrangement events. In the following, we will 
neglect the thermal random force, which would add a Debye- Waller short time 
diffusive contribution to the dynamical structure factor, and set ff = 0. However 
the presence of ff has an indirect crucial effect since the local rearrangements are 
probably thermally activated. 

The slow rearrangement regime 

A first step is to calculate the Fourier transform of the time derivative of the 
displacement field u(f,t) created by a single dipole located at fj, in direction fij, 
that we denote v(k,t\fj,fij,tj). One finds: 



The displacement difference between t and t + r can therefore be expressed as: 



C(q, t) = (exp [iq- (u(f, to + t)- u(r, to))]) 



(56) 



v(k, t\fj,fij,tj) 





(58) 
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The analysis of this expression reveals that there are several cases to consider for 
the relative position of the time tj when the j th event takes place and the other 
relevant times. We refer to |23) for further details. In the course of the analysis, 
one finds that a new, g-dependent time scale T q appears, defined as: 



D9 
( 

qvo 



(59) 



such that, depending on the ratio t/r q , the dominant contribution to the decay of 
C(q,t) comes from different regions of the t',tj plane. In the slow regime where 
the rearrangement time 9 is larger than the lag t, the final result reads: 



C(q, t) = exp 



-A'p9(D9) 



3/2 



3/2' 



(t < T q ), 



(60) 



which has the compressed exponential form found in [22], in particular, it indeed 
only depends on {qt) 3 / 2 . (The numerical value of the prefactor A 1 can be computed 
exactly to be 1.67996...) This result is however only valid in the short time regime 
t <C 9 and t <C T q . Note that the combination p = p9(D9) 3 ^ 2 is adimensional and 
represents the average number of events taking place within a time interval 9 and 
within a diffusion volume (D9) 3 ^ 2 . 



Other regimes 

Other regimes can be analyzed similarly (see [23]), and lead to different scaling 
behaviours for C(q, t). A summary of the results in terms of H(g, t) = — log C(q, t) 
is as follows. The two physical cases depend on the relative position of r q and 
T q = qvo/D = 9 2 /r q . For D9 <C qvo one finds T q <C -C f q and: 

Z-iqtf' 2 (t<r,); H ~ {qtf^ (r,«k<T 9 ); E ~ q 3 ^ 2 t (t»f g ). 

(61) 

Note the appearance of a second compressed exponential regime with (3 = 5/4. 
For D9 3> qvo, this (qt) 5 ^ regime is squeezed out and the results are simply: 

-~~{qtf' 2 (t«C0); E~q 3 ' 2 t (t » 9). (62) 

where we have skipped all prefactors. The late time, purely exponential regime 
in q 3 / 2 t in fact reads: pv^ 2 q 3 ^ 2 t, which can be simply understood. The factor pt 
simply counts the average number of events per unit volume between to and to + 1 } 
and q 3 ! 2 reflects the fact that the distribution of local displacements u decays as 
u~ 5 / 2 and has a diverging variance. For a distribution with a finite variance, one 
would obtain the usual q 2 dependence. 
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The mechanism leading to a compressed exponential at short times is the fact 
that the motion of particles is not diffusive in this regime, but rather ballistic, due 
to a local drift of the structure imposed by slow events taking place elsewhere. 
Clearly, if was very small (corresponding to instantaneous rearrangements), one 
would lose the compressed exponential behaviour altogether and observe simple 
exponential decay with an anomalous q 3 ^ 2 scaling, exactly as for Levy flights (see 
Eq. (|45p ). The fact that the motion of particles is coherent over a rather long time 
scale 9 is crucial to observe the small time slow decorrelation 1 — C(t) ~ t 3//2 . As 
recently noted |24j . this ballistic-like motion might be a generic property of glassy 
system, where "crossing a barrier" means a slow, collective rearrangement of many 
particles in an orderly manner. There might therefore be deep inter-relations be- 
tween compressed exponential relaxation, cooperative effects and dynamical het- 
erogeneities (for recent reviews, see [251 [26] an d refs. therein). 

4 Conclusion 

We have tried to give an overview of the physical mechanisms leading to anomalous 
relaxation, and the relation of this phenomenon with anomalous diffusion and 
transport. Whereas in some cases these two notions are indeed deeply related, this 
needs not to be the case. We have shown in particular that stretched exponential 
relaxation can be found in trap models where subdiffusion occurs, but only as 
a short time approximation, and in a non stationary (aging) situation. At long 
times, the relaxation becomes a power-law in time. We have pointed out that 
stretched exponential relaxation, while being slow at long times, in fact reveals 
the existence of fast relaxation processes at short times. We have emphasized the 
fact that any distribution of relaxation times will lead to a relaxation function 
resembling, at least within some time interval, a stretched exponential; therefore 
one should be very careful about the physical consequences one infers from a 
stretched exponential fit of experimental (or numerical) data points, in particular 
the physical meaning of the stretching exponent (3. 

Models leading to compressed exponentials are much more scarce. In partic- 
ular, the dynamical structure factor of particles undergoing Levy flight motion is 
a simple exponential in time. Similarly, superdiffusion in random advection flows 
lead, at least within some approximations, to a simple exponential in time; it 
would be interesting to know under what conditions this is an exact result. Still, 
compressed exponentials have now been observed in a variety of systems, from soft 
glassy materials to granular packs. We have discussed a generic model for com- 
pressed exponentials where slow rearrangement events occur randomly in space 
and create long ranged elastics strains, leading to f3 = 3/2. 
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